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Abstract 

We calculate the fc-point generating function of the correlated Jacobi ensemble using 
supersymmetric methods. We use the result for complex matrices for k = 1 to derive a 
closed-form expression for eigenvalue density. For real matrices we obtain the density in 
terms of a twofold integral that we evaluate numerically. For both expressions we find 
agreement when comparing with Monte Carlo simulations. Relations between these 
quantities for the Jacobi and the Cauchy-Lorentz ensemble are derived. 


1 Introduction 

The Jacobi ensemble, like the Wishart ensemble, has its roots in the held of multivariate 
statistics. To quantify the empirically estimated canonical correlation coefficients between 
two sets of time series, they are compared to a null hypothesis, i.e. to a Gaussian distribution. 
Assuming for both sets Gaussian statistics with a non-trivial correlation structure, the null 
hypothesis becomes the Jacobi model Besides the canonical correlation analysis it 

applies also to other aspects of high-dimensional inference such as analysis of variances, 
regression analysis and the test of equality of covariance matrices m- Gonsequently, they 
are also referred to as MANOVA (Multivariate ANalysis Of VAriance) ensembles. 

Along with the Gaussian (Wigner) and the Laguerre (uncorrelated-Wishart) ensembles, 
the Jacobi ensemble constitutes the family of classical random matrix ensembles. Just as 
the eigenvalue statistics for Wigner and Wishart ensembles respectively involve Hermite 
polynomials and associated Laguerre polynomials, the eigenvalue statistics of Jacobi en¬ 
sembles involve Jacobi polynomials. Together they complete the random matrix ensemble 
picture in connection with the theory of classical orthogonal, and skew-orthogonal polyno¬ 
mials [6]-[^. Glosely related to these classical ensembles is the less known Gauchy-Lorentz 


1 


ensemble 10,11 . It exhibits a Levy tail and therefore hnds important applications in the 


12 


spectral statistics of covariance matrices in hnance 

Besides their crucial role in multivariate statistics and the intimate connection with the 
theory of classical polynomials, Jacobi ensembles hnd interesting applications in the helds 
of quantum transport and optical hbre communication. In the context of quantum trans¬ 
port, they describe the statistics of transmission (or reflection) eigenvalues for disordered 
mesoscopic systems with ideal leads, see 13 and references therein. The knowledge of these 
eigenvalues, in turn, gives access to important observables such as Landauer conductance, 
shot-noise power and Wigner delay time 14 -19 . This connection of Jacobi ensembles with 


the transmission eigenvalues stems from their relationship with the scattering matrices which 
are modelled using Dyson’s circular ensembles 13,14 . More recently Jacobi ensembles have 


been used in the ergodic capacity analysis of multiple-input-multiple-output optical hbre 
communication 20-22 . In this case the Jacobi structure emerges from the channel matrix 


which happens to be a block of a bigger transfer matrix which is unitary in nature. 

In all these applications listed above one can ask about intrinsic correlations of the 
channels or time series which result in correlated random matrix ensembles. Thus one has 
to relax the condition of the empirical correlation matrix being proportional to the identity 
matrix. This correlation can be attributed to different reasons depending on the context. 
For instance, in time-series analysis data are in general mutually correlated, e.g. see 23 


In the context of multiple antenna communication, this correlation arises because of spatial 
correlation between closely spaced antennas 


24 


as it is the case in any cellphone. To handle 

even 


these situations matrix models with correlations between different time series 
with double correlations 


29,30 , have been introduced; see also Refs. 


24 2N 


and references 

therein. In the Wishart case the level density of the model correlation matrix was computed 
by various techniques. We pursue the ideas of Refs. 26,30 where supersymmetry techniques 
were employed. 

In the case of the Jacobi and the Cauchy-Lorentz ensemble not much progress has been 
made in a similar direction as for the Gaussian because of their non-Gaussian form. Nonethe¬ 
less both models can be traced back to a combination of two random matrices both drawn 
from Wishart ensembles, as explained in detail in sections and 3T To each of the two 
Wishart matrices one can associate an empirical covariance matrix. However, we show in 
section that effectively only one empirical correlation matrix is involved. 

In section we briefly review the Jacobi ensemble before showing the relation to the 
Gauchy-Lorentz ensemble in section The latter is more convenient to apply the projection 
formula and the supersymmetry approach. Nevertheless we present in appendix an 
alternative approach which is a direct way from the integral over the Jacobi ensemble to 
an integral over superspace. The resulting supersymmetric integral is more involved than 
the one resulting from the hrst approach because it involves two supermatrices instead of 
only one. However the second approach exhibits an intrinsic symmetry of the model which 
is not obvious in the more compact result presented in section In section we calculate 
the whole eigenvalue statistics of the correlated complex Jacobi and the correlated complex 
Gauchy-Lorentz ensemble. This is possible due to an underlying determinantal point process. 
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Such an integrable structure is not available for the real case discussed in section For this 
case we need the supersymmetric representation and compute the level density with the help 
of the generalized Hubbard-Stratonovich transformation. Some details of this calculation 
are presented in appendix In section we discuss the level density for the real and 
the complex ensemble in the limit of large matrix sizes in a unihed way. Our results are 
summarized in section 0 


2 Correlated Jacobi Ensemble 


This section is devoted to introduce the theoretical background. The correlated Jacobi 
ensemble is a two-matrix model for a Hermitian matrix H which is composed of two inde¬ 
pendently distributed correlated Wishart matrices and BB\ i.e. 


FFt - 

FFt + HHf 


( 1 ) 


The rectangular matrices F and B have the same number of rows p (the number of time 
series) but may have a different number of columns ni and n 2 (the number of time steps), 
respectively. We assume that ni > p and n 2 > p- In a realistic situation for time series 
analysis both dimensions ni and n 2 are indeed larger than p since the number of time 
steps is usually larger than the number of time series. Moreover the matrix entries are 
distributed as Gaussians with variances [Cf]^^ such that the distribution reads 


P{F\Cf) 


exp(-trFFtG^V7) 


( 2 ) 


and likewise for B, where we introduce 7 = 1 for /3 = 2 and 7 = 2 for /3 = 1. The parameter 
/3 is the Dyson index meaning that /3 = 1 corresponds to real matrices whereas (3 = 2 denotes 
complex matrices. The case of quaternion matrices {(3 = 4) is not considered here but can 
be in general worked out in a similar way. 

The empirical correlation matrices Cp and Cb are hxed as in the discussions of 24 30 
Because of the distribution ([^ the random matrices F and B have upon average the 
covariances as the sample, 


same 


- (FFt) = Cf, - (BB^) = Cb 

rii 712 


(3) 


The angular brackets denote the average over the random matrices F and B. 

We note that instead of the matrix model ([^ one can also consider the matrix model 
FF'^(FF1 -|- BB'^)~^. In that case the eigenvalues he in the interval [0,1]. These two matrix 
models are related by a scaling and a shift transformation. We work here with the matrix 
model ([^ and thus relate directly to the Jacobi weight. Because of the involvement of two 
Wishart matrices in the Jacobi model, it is also known as the double Wishart model. 
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The spectral statistics of H can be studied via partition functions which are averages 


over products and ratios of characteristic polynomials, see e.g. 31,32 

I nil det {H - Ka 2 tp"' 




-P{F\Cf)P{B\C^ 


B) 


■' Hfeil det (-ff - Kbllpy 

Those partition functions are related to the /c-point correlation functions as 

k 


R’^ 


■p,l^ 


(x) = (n 


\i=i 


- ti5{H -Xjtp) 


(4) 


(5) 


sj E I n 


d 


Zp]p{xi + zLie,..., Xfc + tike, x [,..., x',,) 


Lj=±l \j 


^ -*■ 2mp dx', 

)=i J 


with xi,... ,Xk G M. We underline that this dehnition of the fc-point correlation function 
comprises self-energy terms, in particular terms proportional to 6 {xi — Xj). Despite those 
terms the dehnition ([^ is very helpful when the joint probability density of the eigenvalues is 
not explicitly known as it is in the case of the correlated real ensembles discussed in section 
In the case when the joint probability density is known one can apply another dehnition 
of the fc-point correlation function, e.g. see (3^, 




vA^) = = 

\j=i 


{p-k)\p 


pi 


Rpj^{x) — self-energy terms. 


( 6 ) 


where Ei,... ,Ek are the eigenvalues of H. This dehnition is employed to the correlated 
complex ensembles studied in section Both dehnitions are normalized in such a way that 
J dxkRp f^{x) = RAAA ^^d / dxkRp isiA = d'he eigenvalue density of R is given 

by the case k = 1, i.e. 


>5/3(x) := ^^tr(5(hf-xlp)^ = R^ix) = R^ix). (7) 

It is computed in sections and 

The partition function Q is normalized as lim^-^-oo ■^p) 3 ^^(«)det^(—Ki)/det^(—^ 2 ) 1- 
Note that in the case Cp = Cb, the partition function and, hence, the spectral statistics of 
H become independent of the empirical covariance matrices such that the joint eigenvalue 
distribution of H is given by the uncorrelated Jacobi ensemble [^|^[^. Indeed we can rescale 
the matrices E —)■ y/CpE and B —)• \/CfB such that the partition function only depends on 

C'eff = cA'^CbCA^^ i.e. 


ryk\ 1^2 


(k) 



riglidet {H-Ka 2 tp) 
njiidet (iJ- KfeiUp) 


P(F|lp)P(5|C'efr). 


( 8 ) 


Thus in the case Cf 7 ^ Cb, the eigenvalue statistics become non-trivial, even in the complex 
case (/3 = 2), see section]^ Thus, it is reasonable to apply other methods than the standard 
Jack or Zonal polynomial approach [2 34,35 . 
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3 Correlated Cauchy-Lorentz Ensemble 


The correlated Jacobi ensemble is also related to the correlated Cauchy-Lorentz ensemble. 
This relation is established in subsection |3.1 With the help of this relation we express 
the partition function in terms of integrals over supermatrices in subsection 3.2 The 


advantage of the supersymmetric integrals is the drastic reduction of integration variables. 
The parameters ni, n 2 , and p only appear as external parameters in those expressions. This 
fact allows an asymptotic study when these parameters are large which is considered in 
section 1 


3.1 Relation between Jacobi and Cauchy-Lorentz 


To see the relation between the Jacobi and the Cauchy-Lorentz ensembles we take advantage 
of the fact that is generically invertible because rii > p. Then we can rewrite F = FLl 
as a product of a p x p square matrix F which can be either real {f3 = 1) or complex {f3 = 2) 
and a random projection fl G 0(ni)/[0(p) x 0(ni — p)] for /3 = 1 and fl G U(ni)/[U(p) x 
U(ni — p)] for (3 = 2. The measure transforms as d[F] = det*'’^^“^^^'^FFld[F](Jp(n) where 
d/i(n) is the Haar measure induced from those on the groups 0(?7,i) or U(ni), respec¬ 
tively. The projection drops out and the integral over it yields a constant. The same 
procedure can be applied for the matrix B yielding & p x p matrix B with the weight 
exp[- tr C'-ffiFFt/^]det(”2-p)/75^trf[ widehatB] 

Since F is generically invertible we can rewrite any average of an observable O of H in 
the following way 


\ \FFi + BBf)}~\ \FFf + BBfJ 


\ \lp +F-^BB^F-^1) / ' 


(9) 


Note that the observable O has to be invariant under the group action of 0(p) and U(p), 
respectively. Due to the invariance of the probability weight F(F|]lp) the average ([^ over 
F, only, depends on the singular values of B such that it is also true 


O 


FFt -FFt\ 
FFt + BBf) 


= O 


1. 


F-^B^BF-F 


t„ + F-^BfBF-F 


( 10 ) 


This identity is reminiscent of the weak-commutation relation proven in 36 . When rescaling 
B —)■ BF^ the integral over F becomes a deformed Gaussian and reads 


j d[F]detF^+^^-py'fFF^ exp 


■ — tr F 

7 


FtC-g'F + 1, 


F^ 


Thus the spectral statistics of H is equal to the statistics of 

^ I, - F'F't 


Ir 


F'F't ’ 


( 11 ) 


( 12 ) 
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where the p x n 2 matrix B' is drawn from the correlated Cauchy-Lorentz distribution 


pCL 

ni+n2 


(5'IC'eff) = 


) = (t] - P + J + ^)/^] \ 

+ ^^2 - P + j + l)/7] / 


xdet“('^^+’"2)/^ + Ces) . 


( 13 ) 


Indeed the factorization B' = SII' in a square matrix B and a projection If' G 0(n2)/[0(p) x 
0(77-2 — p)] for /3 = 1 and If' G U(772)/[U(p) x U(77 2 — p)] for /? = 2 is still possible. Then we 
would have the additional term det02“Ph7^^t the weight (13). 


Note that the eigenvalue statistics of H' is completely determined by the eigenvalue 
statistic of B'. This means when we calculate the fc-point correlation function of B'B'"^ then 
the fc-point correlation function of H' is given by the substitution b = (1 — x)/(l + x) where 
b is an eigenvalue of B'B'"^ and x is an eigenvalue of H'. Therefore we consider the partition 
function 


= / <'1^1 


-.rciidet (B'-B't - 


nhl det (-S'B't - Kjilp) 


^CL 

ni+n2 


(5'ia 


eff J 


(14) 


Fr om this point on everything works analogously to the correlated Wishart ensemble studied 
. When plugging Eq. (12) into Eq. (|8j) we obtain an explicit relation of 


m 

to 


26 


27 


Tie partition function of the Jacobi ensemble which is 

7fcl|fc2/ \ _ / 1 Ufci-teip + ^ 62 )^ rylkl+k2\k2+kl / 1 “ 


-1 


K2 


k2 : 


1 + 7t2 


-1 


fcl 


(15) 


Especially the /c-point correlation function remains effectively unaffected since in this case 
we have ki = k 2 = k and the additional characteristic polynomials cancel. The fc-point 
correlation function with self-energy terms is then 


k 

i=i 


and similar for the dehnition (|^. We have written the differentials to underline the transfor¬ 
mation properties under changes of coordinates. In the particular case of the level density 
we have 

S'p{b)db= (^ti5{B'B'^ -blp)db^ =S'p ^ (1'^) 

Establishing this relation between the correlated Jacobi model and the correlated Cauchy- 
Lorentz ensemble is the hrst main result of our work. 


Wi-ii5{B'B'^-bjtp)dbj 

\j=i 


(16) 


= R' 


p,P 


1 — X 


n 


2dxi 


1 + Xj fJl (1 + Xj) 


2 


) n 

i=i 
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3.2 Projection Formula 

Before we come to supersymmetry let us refer the reader who is not familiar to superanalysis 
and superalgebra to the textbook by Berezin [^. A general introduction to the supersym¬ 
metry method in random matrix theory is given in and references therein. 


To apply the projection formula introduced for chiral ensembles in Ref. 10 we rescale 
the matrix B' —)■ \/C^B'. Then the partition function (ITdj) is 


j d[B'] sdet-'(R'Rd ^ _ ^--1 ® >^)P^^{B'\t^). (18) 

We use the short-hand notation of the superdeterminant which only encodes the products 
and ratios of determinants. Moreover we understand k = diag(Kii,..., Ki 2 , ..., K,k 22 ) as 


a {k\k) X {k\k) diagonal supermatrix. The representation (18) directly reflects the duality 


between the ordinary matrix space and the supermatrix space. 

Let L = signlm/t be the sign of the imaginary increment in the source. Then we can 


apply the projection formula 10,11 to the partition function (18) which reads 


Z'lf\K) = j d[x]sdei-^/\t,®xX^ (19) 

X sdet"^/^(xx^ -f 

If7fci|7fc2) 

= K2 det^^-'^iC'eff J dfi{U) sdet”^/"^ 17sdef^/^(]lp 0 f/ - ® {Ln O 1 ^,)) 

X sdet“^/'^((L ^1^)U l^ki\'yk2) 

= A3 det^^“^^Cefr y (i[cr]/„2((T) sdef^ (ip ® a — (g) (k (g) 1^)) 

X sdet"^/'^ {a + 17^117^) 

with = rii + n 2 — p + ki — k 2 . The normalization constants are 

= j d[x\sdet~>^^'^{xx^+ ^k^\'fk2), 

dp{U) sdet“'^/^((L 0 t^)U + 17 ^ 117 ^ 2 ) sdet"^/^ A, (20) 


K2^ = 


= 


d[a]In 2 {o') sdet (a + 17 ^ 117 ^) • 


The rectangular supermatrix y has dimension ( 7 ^ 117 ^ 2 ) x ^2 and it satishes the symmetry 

X* = diag(]l.^fc2,r2 0 {x^T = X^diag(l^fci ,-^2 0 Ife) (21) 

in the real case (/3 = 1) where T 2 is the second Pauli matrix, i.e. y consists of a ipfki) x n 2 
real or complex matrix depending on (3 and a ( 7 ^ 2 ) x 77-2 matrix comprising independent 
Grassmann variables. 
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In the second equality of Eq. (19) we have used the superbosonization formula 40 -42 


We use this representation in sections and In the present situation the matrix 


U = 




FF 


( 22 ) 


satishes the following symmetries: 


• [L ® 1^)Ubb = [{L ® ]I^)t/BB]^ positive dehnite for both /3 = 1, 2 and Ubb = 
for /d = 1, 


ul^ = m 


"FF ~ "-"FF unitary for both /3 = 1, 2 and Uyf = T 2 UppT 2 self-dual only for /d = 1, 


• h^BF = f^FB both /3 = 1 , 2 and ffjg = UbfT 2 only for f3 = 1. 

The set of supermatrices U is a. particular case of the co-set Gl{k\k)/\J{k — r,r\k) for /3 = 2 
and UOSP(2fc — 2r, 2r|2/c)/U(2fc — 2r, 2r|2fc) for /? = 1 (r is the number of minus signs in L), 
see 1^. Note that we have chosen the non-compact group symmetries due to the non-trivial 
signs L. 

Alternative to the superbosonization formula one can also choose the generalized Hubbard- 
Stratonovich transformation l42||4S|45], see third equality of Eq. (19), which we employ in 


section The superbosonization formula and the generalized Hubbard-Stratonovich trans¬ 
formation are equivalent [^. The hrst can be understood as the contour representation 
of the latter which is some kind of a high dimensional residue theorem. The distribu¬ 
tion /n 2 (u) in the generalized Hubbard-Stratonovich transformation is the supersymmetric 
Ingham-Siegel integral 


42(u) = / d[p]sdet ’""/^(p-*el^fci| 7 fe 2 )exp [«str(p-?el^fcp..yfc 2 )cr] 




with e > 0. It encodes derivatives of Dirac delta functions in the fermion-fermion block, ctff, 
and the positivity condition in the boson-boson block, ctbb- Therefore the supermatrices 
a and p are elements of another realizations of the cosets G\{k\k)/'[J{k\k) for (3 = 2 and 
UOSP(2/c|2fc)/U(2fc|2/c) for /3 = 1. They are of the form 


a = 

O'BB 

‘^FB 

and p = 

Pbb 

Pfb 


CTFB 

%( 7 FF 


Pfb 

*Pff 


(24) 


where apB and ppB comprise independent Grassmann variables, only, and satisfy Upg = 
'i>T 2 (Jfb and ppg = tT 2 PFB ioT (3 = 1. The blocks (Tbb, upp, pbb, and ppp are Hermitian. 
Additionally the submatrices ctbb and pbb are real symmetric while upp and ppp are self¬ 
dual for /3 = 1. 


4 Eigenvalue Spectrum for [3 = 2 


In the case of complex matrices {(3 = 2) we hrst calculate the joint probability density before 
applying supersymmetry. This simplihes the whole calculation a lot. In particular we show 


in subsection 4.1 that the joint probability density follows a determinantal point process. 


The corresponding kernel is expressed in terms of a supermatrix integral in subsection 4.2 















4.1 Determinantal Point Process 


We start from the weight (13). Let b = diag(6i ,... ,bp) > 0 be the eigenvalues of B'B'"' = 
UbU"^ with U G U(p). The computation of the joint probability density of b is possible 
because the integral 


'U(p) 


d/i(f/)def”^"'^2 [UbU^ + Ces) 


(25) 


( .^p(p-D/2(T^ jKni + n2-p-j-l)\ \ det[{bi + Aj)P^^ ^] i < m<p 

^ ^ 11 in. AJb^AJW 


w=o 


(ni + n2 - p - 1)! 


Ap{b)Ap{A) 


is well-known The matrix A = diag(Ai,..., Ap) > 0 are the eigenvalues of the 

empirical matrix Ces- The measure d^{U) is the normalized Haar measure on the group 
U(p) and Ap{b) = Yli<i<j<pibj - h) = det[&r^] i<ij<p is the Vandermonde determinant. 
Thus the joint probability density is 


*"(6) = /fjpd det ”'b > 0, 


(26) 


with the normalization constant 


A"jpd 


(ui +n2- p)! 


p-i 

P- ("2 - P +-3 + !)■' 


(27) 


The normalization can be readily checked when applying Andreief’s integration theorem 48 


The fc-point correlation function of the correlated Lorentz-ensemble without self-energy 
terms is 


y poo 

r'1,2 (&i,..., ) = n / (b). 


(28) 


j=k+l ■ 


We employ a modihcation of Andreief’s integration theorem derived in 49 and hnd a de¬ 
terminantal point process 


R'lJb,,...,b,) = 


-l)^(p-/c)!iLjpd 


Ap(A) 

(p — k)\p^ 


det 


0 b'^-Pib^ + 

^i-i {n 2 -p + i- l)!(ni 


{rii + n2 - p)\ 




p! 


det [K{bubm-,A)\ 




( 29 ) 


We emphasize that the case fc = p is the joint probability density. The kernel is 


A'( 6 ,, 6 „;A) = - 5 ^(-l) 


P —1 


(rii + n2- p)! 


{n 2 -p + i - l)!(ni - i)\ det(Aj]lp_i - A^j) ^ 
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— ni—722 —1 


with 




4(A) = 5: 

i<ji<...<ifc<p 


Aji... 


( 30 ) 


(31) 


the elementary polynomial of degree k with p arguments and A^j the set of (p—1) eigenvalues 
A where Aj is excluded. Deriving the determinantal point process (29) for the complex case 
is our second main result we have been aiming at. 


The result (29) immediately yields the fc-point correlation function of the correlated 


Jacobi ensemble via the explicit relation 




n 




,/k 


j=i 


X 


(32) 


The prefactor is the Jacobian from the change of coordinates b = (1 —x)/(l + x), cf. Eq. (16). 
In particular the level density ([^ reads 


(ui + n2 - p)\ 




V (77-2 - p + i - l)!(ui - z)! det(Aj]lp_i - A^j) ^ 




X 


(1 + xY Vi + ^ 


fl-x 


+ Aj 

1 + X ^ 


p—ni—n2 — l 


(33) 


This closed-form expression for the level density is plotted in Fig. [T] where we compare the 
result with a Monte-Carlo simulation. The structure of Eq. (33) is reminiscent to the result 


found for the correlated complex Wishart ensemble 24| 29 
Despite the fact that the kernel 


is very explicit it is unsuitable for studying its 
asymptotics, see section This is the reason why we want to rewrite this expression in the 
next subsection. 


4.2 Supersymmetry and the Kernel 


Let us point out that the kernel (30) is independent of k. This is quite convenient since we 


can express the kernel as an integral which is very similar to the one we started from, see 


Eqs. (26) and (28). Thus the kernel is 

roc 


K{xi,X 2 ]A) = iLjpdderiATT f 

.■1 Jo 




' Ap(A) 


(X2 + 


X. 


n2—p 


l<i<p-l 

l<t<P 


(34) 

The normalization constant is fixed via the asymptotics of the kernel in the variables xi and 
X2. 
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Figure 1: Comparison of the analytical expression (33) (solid curve) with Monte Carlo 


simulations (histogram). The random matrices have the sizes p = 3,ni = 5 ,^2 = 7 and 
the empirical eigenvalues are A = diag(l/3, 2,9/2). We have generated 50 000 complex 
correlated Jacobi matrices. 


Note that we now integrate only over p — 1 variables instead of p. The missing integral 
can be introduced by a Dirac delta function 6{bp — X 2 ) = lim£_ 5 .o Im l/[ 7 r( 5 p — X 2 — le)]. The 
symmetrization in all 61 ,..., bp reads 


K{xi,X2;A) = — 

Tip 


1 f 

limlm-Wpddet'^^A IT / 

X2 - Xi Jo 


db. 


Ap( 6 ) det (6 — xi'tp) 
Ap(A) det (6 - {x 2 + ie)tp) 


X det [ 6/^ ^{bi + kjY "2 1 j 


(35) 


I 

What did we gain from rewriting the kernel? We can now identify the integral on the 


right hand side with a partition function (14), i.e. 


K[xi,X 2 ]A) = —limlm 


Z'^^l{x2Aie,Xi). 


Tip e-s-0 X2 — Xi 


(36) 


This allows us to apply the result (|19|) of the projection formula, 

1 


K{xi,X 2 ]A) = lim 


_ - [ da(U)sdet Vlp ® f/— A ^ ® di&g{x 2 + iLe,xi)) 

£-^o 2tiip X2- Xi J 

Li — i i 


X sdetP“”i“”^(LD + li|i) sdet"^ U. 


(37) 
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The explicit parametrization of U is 


U = 


Le^ r]* 


V 




(38) 


with r] and 7]* two independent Grassmann variables and "d G M and (p G [0,27r]. The Haar 
measnre is dn{U) = Le'^~^^'^d'&d(pdr]*dr]/(2'K) where we choose the convention J rjdri = 1. The 
expression (37) is exactly what we are aiming at. The integral is over a small hxed nnmber 


of variables and the “large” dimensions rii, n 2 and p appear as external parameters. Thus 


the result (37) invites for a saddle point approximation. This is the main idea behind the 


supersymmetry method 38,39 


The diagonal elements of the determinant (29) are the level density at positions xi,... ,Xk 
whereas the kernel becomes 


Si(b) = 


lim 

e^o Amp J 


dp{U) str(lp ® 7/ — (6 + iLe)A 


-1 


11 


i|iJ 




11 


0 


p 

0 -1. 


X sdet {lp®U — {b + tLe)A ^ (g) Him ) sdet^ 


0 — 711—712 


{LU + Im) sdet^^ U 


(39) 


due to L’Hospital’s rule. 


5 Eigenvalue Density in the Real Ensemble 

In the case of the real correlated Jacobi ensemble the computation of an arbitrary fc-point 


correlation function is highly non-trivial since the corresponding group integral (25) is un 


known. Hence we concentrate on the calculation of the level density ([^. Already computing 
the level density with the help of the projection formula is quite involved. The dimension 
of the supermatrix model is (2|2) x (2|2) and thus twice as large as for the kernel of the 
complex matrices. 

We apply the generalized Hubbard-Stratonovich transformation |4^|44l|45] for k = 1, 


third equality of Eq. (19), which explicitly reads in this case 

d[a]In 2 {o') sdet“^/^ (a + II212) sdef^^^ (A“^ (g) k 


7/^ 




(dvr) 


tr 


a 


). (40) 


The exponent is p = ni + n 2 —p and is the supersymmetric Ingham-Siegel integral (23) 

for ki = k 2 = 1. Furthermore we have k = diag(6 ±te,b ± te, bi, bi) and 


a = 


(^1 C 

-(i ZCT2112 


c = i; yi '^1 = 1; )). ( 41 ) 

The normalization has been calculated by choosing bi = biLze where Cauchy-like integral the¬ 


orems 


50|-(55 apply. In the parametrization of a, the variables a, b, c, a 2 are real commuting 

(42) 


whereas a, fd are complex Grassmann variables. The measure is 

d[a] = dadbdcda 2 dada* djddfd*. 
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Let us first consider the diagonal blocks. Since the boson-boson as well as the fermion- 
fermion block of K are proportional to I2, we can diagonalize the boson-boson block of 
a without any problems. Accordingly, we write cxi = OrO^, where r = diag(ri,r2) and 
O G 0(2). This change of coordinates yields a decomposition of the differential 


d[ai] = |ri — r2| dridr2d^{0). 


(43) 


Because of the structure of the integrand (40), the integral over 0(2) factorizes and yields a 
factor of 7r/2. 

In the next step, we expand the superdeterminants of cr-|- II212 in the Grassmann variables 
a, a, (3 and /3*, 


sdet (a + I212) = 7--W9 


X 1 + 


[X 


-OL (y. “1“ 




-(3*13 + 


/i(/i - 1) 


(1*(T2)(1ri) {I + ia2){l + r2) {I + ia2Y{l + ri){l + r2) 

det(fi;2A“^ — i(T2^p) 


(44) 

a*al3*l3 


and 

sdet“^/^ (A“^ ® k — Up ® (t) 

/ p 

E 


X 


iydet(KiA“^ — rilp) det(/s;iA“i — r23Lp) 

1 .A 1 


^ (fi:2A^- -ta2)ifXiAj - n) 

1 


a*a + 


^iK2Aj -ia2){niA- - r2) 

a*a(3*(3\ 


-yf3 


('‘2 V “ - *<22)(k 2A, ‘ - !Cr2)(KiAj ‘ - ri)(K,A, ‘ - r^) 


(45) 


The expansion in the Grassmann variables of the supersymmetric Ingham-Siegel integral for 
(3 = 1 was done in Ref. 26 and is 

dvr 


In^ir) = 


(n2-2)! 
X I 1 — z 


0(ri)0(r2)(rir2)^”^ LA | ^ 
(3*(3\ d 


d 


n2—2 


a*a 


n 


+ 


d(T2 
a*af3*f3 d'^ 


r2 J da2 


rir2 dao 




(46) 


where 0 is the Heaviside step function. 

When performing the integrals over the Grassmann variables we only keep the leading or¬ 
der terms in the four Grassmann variables, in particular those terms proportional to a*a(3*(3. 
Moreover we evaluate the Dirac delta function and apply the derivative in K 2 . In the end we 
set K 2 = txi = K. Then the partition function is 


dZ'^ 1 

= g [KA - 1 ) 0 ^ 2 - 2 . 74-2 - 2 /i(rZ 2 “ (k; A) 


xAr°^_^p_+i^_i(fi;;A) +n2(n2 - 1)CP^^^{k; A)] 
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+ z5Z - (rzs - _i(k; A; A^)] 


i=i 


+ i E C''-_\„(fc;A^yX_i,„(K;A;A,,Aj) 


(47) 


l<i<j'<p 


with A^j the diagonal p — 1 dimensional snb-matrix of A removing Aj and A^j j the diagonal 
p — 2 dimensional snb-matrix removing Aj and Aj. The fnnctions Ec) are the 

integrals over cr 2 . The derivatives of the Dirac delta fnnctions can be written as contonr 
integrals and as a finite sum 


= A / ' + e-')* ^(nEj' - 

C 

= E(-i) 




b\ j 


i=i 


r (48) 

(c — a — j)!(a-I -0 — c-1-jj! 


The resulting sum involves the elementary polynomials (31) which already appeared for the 
correlated complex ensembles. The function 




dridr2\ri — r2\ 


ya-l+di)/2 (a-l+d2)/2 


nf=iV«A; -riJKAi^-r2 


.-1 




(49) 


'K nr= 

appears for any term in the expansion of the Grassmann variables as long as we take the 
hrst term in the superdeterminant (45). The second and third term in the expansion (45) 
yields the two-fold integrals 


A; Aj) — 


dridr2\ri — r2\ 


al2 {a-l+d)l2 

'1 '2 


rtp 


-da/'^A 


-1 _ kA- ^ - ri (1 ri)G2(i + r^Yb+i+d )/2 
• 2 '' 


(50) 


and 


<6(«:;A;Aj,Aj) = 


dridr2\ri — r2\ 


niLiV'^A; -tiJkA 


.-1 


r2 


(kA. -ri)(KA. -rs) 


^a/2 a/2 


X 


(l + ri)G2(l + r2)V2’ 


(51) 


respectively. These two integrals are principal value integrals at the non-integrable singular¬ 
ities ri = kA~^ and r 2 = kA~^. We discuss the numerical evaluation of these two integrals 
in appendix [Bj 

The structure of Eq. (47) is very similar to that obtained for the ordinary and the doubly 
correlated Wishart model computed in Refs. 26,30 . The full expression can be separated 
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into three parts; a part with only square root singularities in ri and r2, a part including an 
additional 3/2-singularity in ri, and a part with two 3/2-singularities, one for ri and one for 


r2. The latter two terms in Eq. (47) correspond to the single and double sum. They have to 


be regularized by Cauchy principal value integrals which is done in appendix 

The £ —)■ 0 limit of the imaginary part is taken from the integrals (|4^, (|50|, and (51), 


only. Thereby we adapt the analysis from Refs. 26,30 and concentrate on the product of 


the square roots. We choose the branch cut of each square root along the negative real line. 
Then for any ?/ G M we have 


lim = 

£^■0 y/y — iLe 


Q{y) + LiQ{-y) 


(52) 


Recall that L is the sign of the imaginary increment. Thus the imaginary part of Eq. (47) 
is only non-vanishing in the limit e —)■ 0 if det(6A“^ — rilp) det(5A“^ — r2llp) has a negative 
real part. We recall that n = b + le. Let us assume that 0 < Ai < A2 < • ■ ■ < Ap < cx). 
Then we can divide the integration domain [0, cxd) into disjoint subsets 


|0,«D) = UVj 


(53) 


i=0 


where 


Vo = [0,6Ap ^) , Vp = (6A1 \ 00 ) , and V* = (ftA^ foA^/•) for i = 1,... ,p - 1. (54) 

This decomposition implies det(6A“^ — rip) = (—1)^| det(6A“^ — rUp)| for r G V;. 

We plug the decomposition (54) into the integrals (|4^, (50), and (51) and keep only 


those terms which yield an imaginary part. Then the double integral becomes a sum of 
decoupled one-fold integrals 

^00 ^00 _ r r 

/ dri / (ir2|ri - r2| -)■ ^ sign(/i - ^2) / dri / dr 2 {ri-r 2 ). (55) 

Jo Jo Jv,. Jvi„ 


^1 +^2 E2No+1 


Moreover, we get an additional sign in the sum from the square roots. Assuming ri G 
and r 2 G V/2, we have 



k=i yj{{h + ie)K^ ^ {{h + ie)K^ ^ - r2) 

(_(q+^2-i-i)/2 



:, li -\-12 & 2 N -|- 1 , 


(56) 


r2\ 


otherwise. 


We also obtain for each additional term 1 /(kA - — r) in the integrals (50) and (51) a sign 


sign(p — I — j + 1/2) if r G V). The term 1/2 in the sign function guarantees that the sign 
is positive if j = p — 
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Summarizing everything the level density © of the correlated Cauchy-Lorentz ensemble 


s'M = E / 


dridr2- 


0<ll,l2<p y yy 
«l+« 2 e 2 No + l '1 *2 


-rsl 


k=l 


( 


X 


/i(ri,r2;&,A) 


E 


f3,i,i{ri,r2;b,A) 


V 


f2,i{r 1, r 2] b, A) 

^ ( 6 A^-^ - ri) ' {bA^^ - r,) {bA-^ - rs) 


\ 


(57) 


where /i(ri, r 2 ; &, A), f 2 ,i{ri,r 2 ;b, A), and / 3 ^/_j(ri, r 2 ; 6 , A) are some functions that can be 
read off from Eq. (47). As already said the integrals with the 3/2 singularities have to be 
regularized which is done in appendix]^ Finally we obtain a hnite sum where each summand 
is a product of the three integrals (7^ Ei,Ec) (see Eq. (48)), A) (see Eq. (84)), 

and gl^ ^{b; A; Ai) (see Eq. (g). 




(_l)('.+h+l)/2sign(;j _ 


(58) 


X I fi det 


+ (U 2 - 1) det 


+2 ^ det 
i=i 


0<h,h<P 
/i+/ 2 £ 2 No +1 

0 9n2+l,p+2,hib'^ 9n2-l,p+2,h(b'^ 

{n2 - ( 6 ; A) ^°,+i,^+ 2 ,i 2 (b; A) (b; A) 

. (a - ^)C^,-2,p-2ib; A) gn,-i,^,,ijb; A) <_3,^,)Jfe;A) _ 

0 9n2-l,p,hib'^ 9n2-3,p,h(b-^ 

9n2-l,p,l2(b-^ 9n2-3,p,l2ib-, 

n2CP^,f,{b-, A) ^° 2 +i^^+ 2 ,/ 2 (^;'^) 9n2-l,p+2,l2ib-,^) _ 

0 9n^+i,i,,hib; A; Aj) gl^_^^^ i^{b; A; Aj 

(n2 - l)C'^/_\,^(&; A^j) ^°, + i,^+2,Z2(^; A) 9n2-l,p+2,l2(b'^ A) 

AC:;\,.Ab; A^,) A) <_3,,, A) 

^n 2 +i,/.,q 9^2-1, P,h (b-, A; Ai) 


+2 E ^C-2,.(&;A^7.)det 


1<2<J'<P 


. 9n,+i,f,,i, {b; A; A,) (&; A; A,) 


The integrals g^^c,Ab] A) and gl^c,iib] A; A*) are one-fold integrals over the compact interval V) 
which is numerically more advantageous than the original two-fold integral. 

The level density ([^ of the correlated Jacobi ensemble readily follows from (58) via the 
relation 


S'i(x) 


(1 + x)' 


- S' 

I 


1 — a;\ 
l + x) ' 


(59) 


This third main result is compared to Monte Carlo simulations in Fig. The excellent 
agreement validates our calculation. 
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Figure 2: Comparison between the result (59) (solid curve) and Monte Carlo simulations 


(histogram) for ui = 5, U 2 = 5 and p = 2 with Kj = 1, 4. The sample size consists of 50 000 
real correlated Jacobi matrices. 


We emphasize that the integrals drastically simplify when the spectrum of and thus 
A is doubly degenerate. Then the limit £ —)• 0 in combination with taking the imaginary 
part yields Dirac delta functions such that one of the integrals over ri and r 2 can be exactly 
performed. The remaining integral can be easily performed, too, by employing the idea 
of 49 where the term 1/ det(6A — rtp) can be written as a quotient of a Vandermonde 
determinant and a Cauchy-Vandermonde determinant. The degeneracy of the spectrum is 
not that academic. Such a degeneracy was proposed for the correlated real Wishart ensemble 
to calculate explicit analytical results. It was shown that in the limit of large matrix 


m 


27 


dimensions such an artihcially introduced degeneracy has no effect on the spectral properties 
like the level density and the fc-point correlation function. The open question is if this 
statement carries over to the Cauchy-Lorentz and the Jacobi ensemble. For the level density 
we answer this question in the next section. 


6 Asymptotics of the Level Density 


To obtain the asymptotic behavior of the level density for both, the real and the complex 
correlated Jacobi models, we perform a saddle point approximation of the supersymmetric 
expression (19) for k = 1. Particularly we hrst consider the level density of the correlated 
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Cauchy-Lorentz ensemble 




Tip e-s-0 

K 


^'pA^ + ?£, 6 - J + te) 


j=o 


lim Im / dp{U) exp 


--AU) 

7 


Tip'y e^-O 

The normalization constant K is 

K~^ = [ dp{U)exp 


^ str (6 + te — AjU) ^ 
i=i 


■-m) 

7 


0 

0 0 


(60) 


(61) 


where we nsed an apparent 6 -dependent version which is more convenient in the saddle 
point approximation. The 6 -dependence is indeed only apparent becanse of the Canchy-like 


integration theorems in snperspaces 50 55 


The “Lagrangian” in the exponential fnnction is 


C{U) = str 


ni + n2-p 

P 


In (U + — —Inll + - ^In ((6 -|- ie)AA'^'y\'y — f/) ) . (62) 

^ ^ k=i 


This Lagrangian has to be minimized if rii oc n 2 oc p 3> 1, meaning that its first derivative 
has to vanish, 


C'{Uo) = 


rii +n2 -p 
P 


(C/„ + --jup --Y. (‘Ay 1* - C^o)' = 0 , (63) 


k=l 


Since this eqnation is invariant nnder the snpergronp UOSp(2|2) for /? = 1 and U(l|l) for 
(3 = 2 we may diagonalize Uq snch that each of its eigenvalnes satishes 


C'{qo) = 


ni+n2-p 1 


n2 1 1 Aj 


P 


go + 1 


1 V I: 

n h — 


P go P^A~ 


= 0 . 


(64) 


We perform the same asymptotic analysis as in 27,and connt p -|- 2 poles at Qq = 
— 1,0, 6 A“^,..., 6A7^. The asymptotic behavior of C'{qo) at go = ^A^, ■ ■ ■, bA^^ implies 
that p — 1 of the p -|- 1 solntions of the saddle point eqnation (64) are real. When taking the 
imaginary part in Eq. (60) we recognize that the real solntions do not contribnte. Thns we 
are looking for the complex conjngate pair which solves Eq. (64). 

To hnd a closed form of the saddle point solntion we nnderline that all poles and all 
zero points, apart from the complex conjngate pair, of >C'(go) lie on the real line. Thns an 
integral along an appropriate contonr Cn of its logarithmic derivative yields the solntion via 
Canchy’s integration theorem, i.e. 


go = lim 

R^oo 


'Cn 


9 In C' 

'' dz 


dz 
^ 2tii 


(65) 
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The contour is chosen as follows 


6*^= + [0,7r]}u|r + ^|re 


( 66 ) 


and we integrate clockwise. In this way we obtain either the unique complex solution in the 
complex upper half-plane or zero. The latter implies = 0, too. 

The fermion-fermion block encircles the pole at the origin but the other poles of the 
Lagrangian are purely zero points of the integrand. Hence this contour can be deformed to 
go through both complex solutions go and q^. However the boson-boson block is hindered 
to reach both solutions since the poles of C are apart from the one at origin poles of the 
bosonic integrand. The sign in front of the imaginary increment te dictates which pole can 
be reached and which not. By a contour deformation we can include the saddle point in 


the upper or lower complex half-plane. In Ref. 57 it was shown that the leading order in p 


is given by those saddle points where the boson-boson block and the fermion-fermion block 
share the same saddle point. Thus, in the vicinity of the saddle point we have 


U = (Re go L^Imgo) + 


6U\ 


(67) 


with 6U a Hermitian ( 7 I 7 ) x ( 7 I 7 ) supermatrix satisfying certain symmetries in the case 


= 1, see 42,45 . We plug this expansion into Eq. (60) and integrate over 6U. We arrive at 




E75- 


AjImgo(&) 


( 68 ) 


p ^5 + ^£-A^.go( 5 ) vrp ^ (6 - AjRego( 6 )) 2 -F (AjImgo( 6 )) 2 ’ 

where we already hxed the normalization by the integration J S'p{h)db = 1. As in Ref. 
we can simplify the expression ( 68 ) with the help of the saddle point equation (64). Tnen 


27 


we hnd our fourth main result for the correlated Cauchy-Lorentz ensemble 
ni+n 2 -pl^ qo{b) ni + n 2 -pl Imgo( 6 ) 


SUb) 


Tip 


pi 

— -Im- ,,, 
b qo{b) + 1 


Tip b {Reqo{b) + 1)'^ + {lmqo{b)y 


and for the correlated Jacobi ensemble 
2(ni -f n2 - p) 1 


S^ix) 


Imgo([l -x]/[l -Fx]) 


1 - x2 (Rego([l - x]/[l x]) 1 )^ (Imgo([l - x]/[l x]))^ 


Tip 


(69) 


(70) 


This solution is slightly more involved compared to the correlated Wishart ensemble 27 


Nevertheless it is a closed form in combination with the saddlepoint solution (65). 


We compared the result (70) with Monte Carlo simulations in Fig. We hnd a perfect 
agreement for a matrix size p = 32. Due to the weaker level repulsion of real matrices 
compared to complex ones the agreement is in the bulk better for /? = 1 than for /3 = 2 while 
it is worse in the tails at the edges. Indeed the strong oscillations for the complex ensembles 
is a direct result of the stronger level repulsion which is quadratic. Even the positions of 
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Figure 

value 


3: Comparison of the analytical 

density (solid black curve) with 


expression (70) for the limiting eigen- 
numerical simulations with correlated 
real Jacobi matrices (blue histogram) and correlated complex Jacobi matrices 
(red histogram). Each ensemble consists of 10® matrices with the parameters 
p = 32, ni = 71,^2 = 68 . The empirical eigenvalues were independently and ran¬ 
domly drawn from a Gaussian, in the present case the diagonal correlation matrix is 
A = diag(294.845, 34.679, 30.311,11.612,10.733, 9.468, 8.232, 5.307,4.144, 2.443, 2.429, 2.218, 
2.083,1.986,1.406,1.382,1.102,1.001, 0.889, 0.707, 0.693, 0.684, 0.665, 0.63, 0.594,0.591,0.576, 
0.574,0.562,0.467,0.463,0.455). The largest eigenvalues are three outliers which can be 
seen in the vicinity of a; = — 1 , see the inset. 


the outliers are astoundingly good predicted and the shape of their distributions moderately 
approximated. 

We want to emphasize that we obtain the well-known results for the uncorrelated 
case. Then the distribution of the Jacobi ensemble has either square root zeros or square 
root singularities at the boundaries of the support depending on if the edges detach from the 
generic bounds x = ±1 or not. This can be easily checked by our results. For the Cauchy- 
Lorentz ensemble we obtain a Levy-tail with an algebraic decay of 6 “®/^ if the corresponding 
Jacobi ensemble does not detach from the lower bound x = —1 which was also found in 12 . 

Finally, let us come back to the case of a degeneracy of and, hence, A. When taking 
p —)■ Ip, ni —)■ Irii, n 2 —)■ ln 2 , and A —)■ A(8)]l;, we immediately notice that the number of copies 
I G N drops out in the saddle point equation (64). Thus the saddlepoint is independent of 1. 


Also in the level densities (|69|) and (|70|) the number I does not appear. Hence the asymptotic 

Jistics remain unaffected, too, as it is the 
For the complex case this can be readily 


result is independent. We expect that the local statistics remain unaffected, too, as it is the 
case for the correlated Wishart ensemble, see 


checked via the kernel (37). 
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7 Conclusions 


We considered the correlated Jacobi and Canchy-Lorentz ensemble. Onr first resnlt is the 
map of the correlated Jacobi ensemble to the correlated Canchy-Lorentz ensemble. The 
eigenvalne statistics of the hrst is completely determined by the latter and vice versa. This 
is impressive since the hrst ensemble exhibits a spectrnm on a compact interval namely 
[—1,1] while the second one has a Levy tail which drops off in the nncorrelated case with the 
algebraical behavior The nncorrelated Canchy-Lorentz ensemble was tried to apply to 
hnance 


12 becanse of the known heavy tail behavior in the statistics. 


We also derived a snpersymmetric integral for the fc-point generating fnnction of the 
correlated Jacobi ensemble and the correlated Canchy-Lorentz ensemble via the projec¬ 
tion formnla flOlIlll in combination with generalized Hnbbard-Stratonovich transforma- 

The resulting integral over the 


42 44 j:. 


and the superbosonization formula 40 42 


26 . This reple¬ 


tion 

supermatrices looks similar to the one of the correlated Wishart ensemble 
sentation is ideal for studying the asymptotic limit (p ~ ni ~ n 2 —)■ C)o) due to its small and 
hxed number of integration variables. For example we calculate a closed expression for the 
limit of the macroscopic level density. The approximation is already very good for moder¬ 
ate matrix size p ~ 30 and generic empirical hxed correlations which is conhrmed by Monte 
Carlo simulations. We underline that even the outliers are predicted by the asymptotic limit. 

The level densities at hnite matrix dimensions were explicitly calculated for both the 
correlated real and complex random matrix ensembles. As in the Wishart case the real en¬ 
sembles are more involved. Nevertheless we could simplify the result to a hnite sum where 
each summand is a product of three integrals and one of the integrals can be performed 
exactly. Both technical properties resemble the results for the Wishart ensemble, see |^ . 
The remaining one-fold integrals can be numerically evaluated. In the case of double de¬ 
generacy of the spectrum of the empirical covariance matrix all integrals can be analytically 
performed. This generic degeneracy was proposed in 27 where it was shown that the asymp¬ 
totic spectral statistics of correlated Wishart ensembles do not diher from the case without 
degeneracy. We also observed that the degeneracy has no inhuence on the asymptotic statis¬ 
tics of the correlated Jacobi ensemble and the correlated Canchy-Lorentz ensemble. Thus 
we propose also to study the case with a doubly degenerated covariance matrix artihcially 
introduced by taking two copies of the covariance matrix in the case of real ensembles due 
to its analytical advantage. 

In the real case we had to restrict our explicit calculation to the level density. Higher 
order correlations are not analytical feasible for hnite matrix dimensions at the moment due 
to the lack of knowledge about certain group integrals. This is not the case for the complex 
ensembles where we took a different approach as for the real matrices. For the correlated 
complex Jacobi ensemble and Canchy-Lorentz ensemble we hrst derived the joint probability 
density of their eigenvalues. This was possible due to an Itzykson-Zuber-Harish-Chandra- 
like group integral derived in 46,47 . The joint probability density satishes a determinantal 


point process and the corresponding kernel resembles the result of the correlated Wishart 
ensemble 24,29 . Nonetheless we derived also for this kernel a supersymmetric integral 


which is much more suitable to study the behavior at large matrix sizes. 
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Our analysis can be extended into various directions: First of all, the supersymmetric 
expression for the /c-point correlation function is a perfect starting point to derive also for k > 
1 closed expressions. The distributions of the largest and the smallest eigenvalues are other 
important quantities which were already studied for the correlated Wishart ensemble 25 


28,59 . Also the case of double correlations, as discussed in Ref. 30 for Wishart ensembles. 


could be considered. Another generalization could be the investigation of other correlated 
heavy tailed ensembles instead of the Cauchy-Lorentz one. Product matrices, see 58 for a 


recent review, yield a new and analytical feasible approach to such heavy tailed ensembles. 
The combination of the approach applied in the present work with product matrices and the 
projection formula 10,11 may provide a unique and ideal tool to study the macroscopic as 
well as the local spectral statistics of Levy tailed ensembles. 
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A Supersymmetric two-Matrix Model 


Let us derive another supersymmetric integral which consists of two supermatrices entering 
in a symmetric way. This integral explicitly shows the symmetry of the correlated Jacobi 
ensemble under rii <H- 77-2, A —)■ A“^, and x —)■ —x. This symmetry is not immediate in the 
expression (19) where we have to substitute b = (1 — x)/(l + x) since it is given for the 
correlated Cauchy ensemble. However the supersymmetric integral (19) is certainly simpler 
to compute than the one we present in this section since we have to deal with only one 
supermatrix in Eq. (19). 

We start from Eq. (|^. To apply the same approach as in 
arguments of the characteristic polynomials in the FF^ and BB'. 


we have to linearize the 
This can be achieved by 


multiplying the matrices in the determinants from the right with (FRl + RRl) yielding 


det ((FFt _ 5 Rt)(FFt + _ A + det {FF^Ka 2 - BB^) 

det((FFt-RRt)(FFt + RRt)-i-K^Ap) A + det - RRt) ’ ^ > 

with 'k = (1 — fi:)/(l + k). In the next step we plug Eq. ( [7T| into Eq. (|^ and express the 
determinants as a Gaussian integral over a rectangular supermatrix 


^ ;3 = 1. (72) 

<3 = 2, (73) 
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of dimension p x ( 7 ^ 17 ^), i.e. 

n = /'*'"‘1 “P (' +,trBB>AA^). (74) 

Here we assume for simplicity that the imaginary parts of ki = diag(Kii,..., Kki) are on the 
complex upper half-plane. The source matrix is 


j = diag (Kii, . . . , Kkl,Kl 2 , ■ ■ ■ , Kfc 2 ) <8 


(75) 


We substitute the integral (74) into the generating function (|^ and exchange the F and B 
with the A integral. The resulting F and B integrals are Gaussian and yield 


j d[F] j d[B]P{F\lp)P{B\A)exp {ttiFF^AjA^ + ttiBB^AA^) 
= def'^i/^ {Ip - zAjA^) def”^/^ {ip - zAAA^) . 

Then the generating function becomes 

Z^^{k) oc sdet“^ (Ifc + k) f (i[H]det“"'^'^'’' (ip — iA}A^^ det”"'^^'^ (ip — lAAA^^ 


(76) 


(77) 


The normalization constant is independent of k and A. The next step is known as the duality 


between ordinary and superspace. Due to the invariance of the integrand in Eq. (77) under 


A ^ UA for an arbitrary U G 0(p) for /? = 1 and U G U(p) for /3 = 2, the integrand 
only depends on the invariants tr {AA^^ for m G N. These invariants are equal to the 
superinvariants tr {AA\^ = str see Employing this duality in the 

generating function (77), we arrive at 


Zpi^{k) cx sdet ^ (Ifc -|- k) / (i[A]sdet sdet {\^^^^ — *A^AA) . 

(78) 


The main difference of Eq. (78) to most models discussed in the literature so far is that 


this one includes two different products of A and Ah Namely, A^A which arises naturally if 
invariant matrix models are considered and A^AeffA appearing due to a non-trivial correlation 
structure. We cannot replace both products by one supermatrix, but we can apply the gen¬ 


eralized Hubbard-Stratonovich transformation 42,44,45 independently for both products. 


It yields the following supersymmetric two-matrix model 

Zp^^(fi;) cx sdet“^ (Ifc-|- k) J (i[(T](i[p]J„ 2 (p)/„j(cT) exp (—str p — str 

X sdef^A ® a — A~^ ® pj) , 


O' 


(79) 


where the function ImiQ), * = 1,2, is the supersymmetric Ingham-Siegel integral (23). The 
{'jkl'yk) X (yfclyfc) dimensional supermatrices p and a have the same symmetries as the 


supermatrix in the third equality of Eq. (19). 
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We can completely symmetrize the integral in p and a by going back to the sources k ^ k 

_ 1/2 _ 1/2 

and the empirical matrix A —)■ Cp ' CsCp ' . Then we have the hnal result 


^ / '^Mc^[^’]4i(^')42(cr)exp (-strp - str 


(T 


X sdet {Cb ® <7[(1 + k) (g) l^] — Cp® p[(l — k) ® l^]). 


(80) 


This expression is completely invariant under the original symmetry ni o ^ 2 , A —)■ A and 
K ^ —K because the symmetry is achieved by the change p -H- a. 


We again underline that the supermatrix model (79) can be in principal computed by 


expanding the integrand in the Grassmann variables and performing the remaining integrals. 
However we have now two supermatrices such that this calculation can be a highly non-trivial 
task. This is the reason why we use more advanced techniques which involve the relation 


to the correlated Cauchy-Lorentz ensemble, see section 3.1 The supersymmetry result (19) 


can be obtained from Eq. ( |79[ ) by rescaling a ^ pa and then integrating over p which yields 
the superdeterminant sdet+ 't^k\'yk)- 

B Regularizations of the Integrals in Section 


The numerical evaluation of the integrals (|50j) and (jSiJ) suffer by the non-integrable singu¬ 
larities of order 3/2 at the boundaries, in particular they are of the two forms 




bA 


Jl = 


dr- 


/W 


j-i 


bA 


\bAA 


13/2 ’ 


J 2 — 


dr- 


fir) 


■j+i 


bA" 


|6A^- ^ — r 


3/2 


(81) 


for certain real valued functions /(r) without singularities in the interval [bAjA,bAj^]. As 
already said the integrals are taken via Cauchy’s principal value because of the original 
imaginary increment le. Thus we can effectively regularize the integral as follows 


bA. +e 


Ji = limRe 
£—^0 


dr- 


fi.r) 


bA 


({b + ie)Aj ^ - r) 


3/2 


i+i 
bAJ^+e 


= limRe 
£—^0 


. m - HbAp) 

dr 


bA-/i 


[{b + te)Aj ^ - r) 


+ lim Re 


3/2 £—>0 


bAT 




2fibA. 


- 1 ^ 


bA"/, 

j+i 


|6A7^ 


| 3/2 


|6A 


-1 

j 




'^We thank Petr Braun for showing us this technical trick. 


2 /(&A^ 




^JJb~+ie)A~AAp 


bA. +e 


r=bA 


i+i 


(82) 


24 

















and similar for the other integral (then the imaginary part is needed) 

6 AT 


V-i 


J, = 


dr 


fir) - fihk-^) 2f{bAj^) 


6AJ 


I bAj ^ — r 


3/2 




(83) 


The minns sign in Eq. (83) in front of the second term resnlts from taking the imaginary 
part despite it is evalnated at the npper bonndary. The cnt-off of the intervals can be also 
chosen independently of e. The reason is that the other bonndary term of the integration 
by parts vanishes dne to taking the real or imaginary part, respectively. 

We dehne the following two one-fold integrals 

r°‘/‘^dr 


9lcAb;A) = 


and 




sign(p — I — i) 


Ivi (1 + det(6A“^ - rlp)| 


ra/2dr 


(84) 


(85) 


'vi (1 + ry/^lbA- ^ - r|y^|~det(6W-^Ar7f^’ 

I p - i,p - i + I, 


2 (&A-^) 


-l\a/2 


dr 


X 


= 


(1 + 6A,)=/y|6A-‘ - det(6A;‘ - tA"'VOI 


r.a./2 


(1 + rYA^\dei{bA^] - r VOI (1 + bA^ yA^\dei{bA^] - bA -^ VOI. 

l = p-i, 


2{bAr^) 


-l\a/2 


dr 


X 


(1 + 6A,)=/y|6A-‘ - tApiiyi det(6A;‘ - 6A-‘lp_i)| A-.+. " ’’I’'" 

(6ArA“/2 


r,a/2 


(1 + r)'=/2.^| det(6A^| - rVOI (1 + 5A. i)‘=/2.^| det(6A^,^ - 6A. ^ VOI. 

/ = p — i -|- 1 . 

Then we can combine the discnssion abont the splitting of the integral over into disjoint 
sets and the regnlarization of the 3/2-singnlarities. Therefore we explicitly have for the 
imaginary parts of the integrals (49), (1^, and (51) 


( 86 ) 


ilimImW' j j (* + !£; A) = - J]) (-l)<'i+6+i)/2gign(;j _ 

TT e^yJ ? 5 x, ^ ^ 


0<li,l2<P 
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X det 


fi'a+(ii+l,c+(ii + l,Zi 9a+di-l,c+di+l,li 

_ 9a+d2+l,c+d2 + l,l2^^^ 9a+d2 — l,c+d2+l,l2^^^'^'^ 


-Mm lmNl A-Aj) = - Y1 

71 e^O 71 ^^ 

0<h,h<P 


-l)('i+'"+‘)/^sign(!i - I2) 


( 87 ) 


/i+/ 2G2 NoH -1 
4 


limImiV^^(K; A; Aj, Aj) = 


71 £-S>0 


X det 

1 
71 


9l+2,c,hip] A] Aj) gl^^i^{h]A-,Aj) 

fi'a+(i+l,c+d+l,«2(^4) fi'a+d-l,c+d+l,i2(^4) 

0<h,h<P 

/1+/2G2N0 + I 

^ r 9l+2,c,h (^; 9ic,h (^; 

<c.z2(&;A;A,) 


(88) 


These results can be combined with Eq. (47) to hnd the level density S[{b) of the correlated 
Lorentz ensemble, see Eq. (58). 
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